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BaRuOs having five different crystal structures tias been synthesized by varying the pressure while 
sintering. Contrary to the other phases being nonmagnetic, the cubic perovskite phase synthesized 
recentiy shows itinerant ferromagnetic character. We investigated this ferromagnetic BaRuOa using 
first principles calculations. A few van Hove singularities appear around the Fermi energy, causing 
uimsually high magnetovolume effects of AAf/Aa~ 4.3 /is/A as well as the Stoner instability 
[IN{0)~ 1.2]. At the optimized lattice parameter a, the magnetic moment M is 1.01 /is in the 
local spin density approximation. When spin-orbit coupling is included, the topology of some Fermi 
surfaces is altered, and the moment is reduced by 10 % to a value very close to the experimentally 
observed value of ~ 0.8 fin- Our results indicate that the ferromagnetism is induced by the Stoner 
instability, but the combined effects of magnetovolume and spin-orbit coupling determine the net 
moment. In addition, we briefiy discuss the results of the tight-binding Wannier function technique. 

PACS numbers: 71.20.Be, 71.20.Dg, 75.50.Cc, 75.80.-|-q 
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I. INTRODUCTION 

Ruthenates, which possess formally tetravalent Ru 
ions of Ad^, show interesting properties: unconventional 
superconductivity in Sr2Ru04,'lj the metal- insulator 
transition in La2Ru207 and Li2Ru03,[2|, Q metam- 

] agnetism in Sr3Ru207,[3| and complicated magnetic 
behaviors in double pyrochlore i?i?2Ru207 (i?_E=rare 
earth element). |5J] In particular, the Ru-based per- 
ovskites ^RuOs (^=alkaline-earth elements and Pb) 
have been studied extensively for several decades due 
to their atypical and controversial electronic and mag- 
netic behaviors. 'g'-'T?] SrRuOa is a metallic ferromagnet 

] with Tc=160 K and an experimentally observed mo- 
ment of ~1.6 fiB, whereas CaRuOa is a nonmagnetic 
metal, [l^ [l3| PbRuOs seems to exhibit incipient mag- 
netism, but whether this system is a metal or an insu- 
lator at low T is still under debate, [l^ [l^ The crys- 
tal structure of BaRuOa depends strongly on the tem- 
perature and pressure during synthesis, [17| whereas the 

[ other compounds are orthorhombic. Depending on the 
amount of corner- and face-sharing RuOg, the struc- 

I tures of BaRuOs are categorized as the four-, six-, and 
ten-layered hexagonal phases (4iJ, QH, and lOH, re- 
spectively), in addition to the nine-layered rhombohe- 
dral phase (9i?).[18-2^ As discussed both theoretically 
and experimentally. [I2I ITsl [23j the lOH and 9R phases 
seem to be semi-metallic, whereas the rest are nonmag- 
netic metals. Early in 2008, Jin et al. synthesized the 
cubic phase using the technique of high T (~ 1000 °C) 
and high P (- 18 GPa).^17] 

In this research, we will focus on the cubic phase, which 
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shows ferromagnetic behavior with Tc= 60 K.^^ On the 
basis of resistivity measurements, which show a sharp 
peak in the T-derivation of the resistivity at Tc, Jin et 
al. suggested simple metallic characteristics in the high 
T regime, but Fermi liquid behavior in the low T regime. 
The effective magnetic moment is 2.6 /is, close to the 
value for S=l. However, the saturated moment of 0.8 
/is is significantly less than the value corresponding to 
5=1, implying atypical magnetic behaviors. Moreover, 
when pressure is applied, Tc decreases, reaching 50 K at 
P=3 GPa. 

The elastic properties and pressure effects have been 
theoretically investigated to some extent, [l^, HI] but no 
detailed studies of the magnetic and electric properties 
of this cubic phase have been reported yet. Thus, we will 
address the magnetic properties of this system through 
first principles calculations, including the local spin den- 
sity approximation (LSDA), the LSDA plus spin-orbit 
couphng (LSDA+SOC), fixed spin moment (FSM), and 
tight-binding Wannier analysis. 



II. STRUCTURE AND METHOD OF 
CALCULATION 

First, we optimized the lattice parameter a in the cu- 
bic phase. LSDA calculations lead to our optimized a= 
3.938 A. This value is about 1.7% smaller than the ex- 
perimentally observed value at room temperature, [l^ We 
used the optimized lattice parameter in our calculations, 
since the total moment at the optimized a is close to the 
experimental value, as discussed in the next section. 

Our calculations were carried out with LSDA and 
LSDA-I-SOC implemented in the all-electron full- 
potential code WIEN2k. [2^ The basis size was deter- 
mined by RmtKmaa;=7 and APW sphere radii of 2.5 
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FIG. 1: (Color online) Comparison of nonmagnetic total den- 
sities of states (DOSs) per formula unit between cubic and 9i? 
BaRuOa. For the cubic phase, the optimized lattice param- 
eter was used. In the cubic phase, two sharp peaks appear 
around the Fermi energy Ef, implying magnetic instability. 
The vertical dashed line indicates Ef, which is set to zero. 



for Ba, 2.0 for Ru, and 1.7 for O. The Brillouin zone 
was sampled with 14x14x14 A:-mesh. For both the 
LSDA-I-SOC and FSM calculations, (27j a much denser k- 
mesh of 20x20x20 was used. The lattice parameter was 
relaxed until the forces were smaller than 2 mRy/a.u. 
Furthermore, the all-electron full-potential local-orbital 
code FPLO was used for the tight-binding Wannier func- 
tion analysis. [2^ 

III. RESULTS 

Figure [T] shows the overlapped total densities of states 
(DOSs) between the nonmagnetic 9i? and cubic phases 
between -1 and 1 eV. The DOS of the cubic phase is five 
times larger than that of 9i? due to sharp peaks at Ep 
and 30 meV, suggesting strong magnetic instability. 

Using LSDA, we calculated the change in the total 
magnetic moment while varying a from the experimental 
value to a value slightly smaller than the optimized one. 
As shown in the inset of Fig. [2l the magnetic moment 
decreases linearly with a slope of AM/Aaw 4.3 /is/A, 
indicating strong magnetovolume effects. This value is 
much larger than that of ~1 /xg/A in the Invar systems 
^Fe2 (^=Y, Zr, Hf, and Lu).[2| 

Fixed spin moment studies. To investigate these be- 
haviors, FSM calculations were conducted. Plots of en- 
ergy versus total magnetic moment M for varying values 
of the lattice parameter a are shown in Fig. [51 The dif- 
ference in energy between the AI—0 and ferromagnetic 
(FM) ground states decreases monotonically from ~ 32 
meV at a=4.01 A to 17 meV at 3.93 A. The curves are 
fitted well by E{M) - E{0) w -aAP + jSM"" with con- 
stants a and /3 in low M regime. For the bare suscepti- 
bility xo = 2/i^A^(0), the Stoner-enhanced susceptibility 



FIG. 2: (Color online) FSM calculations in the cubic phase, 
while changing the lattice parameter a (in units of A). E{Q) 
denotes the energy of the nonmagnetic state for a. The arrow 
roughly connects the minimum energy states. Inset: Change 
in total magnetic moment M with respect to a in the FM 
calculations. The arrow indicates the value at the optimized 
a . A discontinuity occurs at a ~3.91 A, since Ef pins a van 
Hove singularity in the spin up channel (see below). 
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FIG. 3: (Color online) Change in Stoner parameter / (top), 
DOS N{0) per single spin at Ef (middle), and IN{0) (bot- 
tom) versus lattice parameter in the cubic BaRuOa. Roughly, 
IN{0) decreases linearly as the lattice parameter decreases. 



is given by x = xo/[l — ^i^)!]: where A^(0) is the single- 
spin DOS at Ep. In low M regime, a = ^X~^- The 
top panel of Fig. |3] displays the obtained Stoner param- 
eters / w 0.70(±0.02) eV for each a. As a changes, the 
DOS at Ep fluctuates substantially owing to van Hove 
singularities (vHSs) near Ep (see below), as given in the 
middle panel of Fig. [S] IN{0) decreases monotonically 
as a decreases, which is consistent with the substantial 
magnetovolume effects, but for all the lattice parameters 
studied here these values exceed the Stoner instability 
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FIG. 4: (Color online) Total and atom-projected DOSs of the 
FM cubic BaRuOs in the full range, including all O 2p and 
Ru 4d states. In the spin up channel, peaks appear at -0.2 
and -0.1 eV. DOS N(0) at Ef is 3.34 states per eV for both 
spins. Upper inset: Expanded view of total DOSs in LSDA 
and LSDA+SOC between -0.5 and 0.5 eV. Lower inset: Ru 



t2g and eg orbital-projected DOSs. 




FIG. 5: (Color onhne) Blowup FM band structure in LSDA 
(top) and LSDA-I-SOC (bottom) near Ef in the cubic 
BaRuOa. In the top panel, the solid and dashed lines rep- 
resent the spin up and down characters, respectively. Around 
the X and M points, the effects of spin-orbit coupling are 
substantial. The horizontal dashed lines indicate Ef- 



criterion IN{0) > 1. 

Electronic structure of FM. Like the FSM results, the 
LSDA results indicate that the FM state is energetically 
favored over the nonmagnetic state by 16.4 meV at the 
optimized lattice parameter. For comparison, from a 
simple Stoner instability, the energy gain of Jm^/4 with 
/ «0.7 eV is 175 meV, much larger than the calculated 
value. 

Now, we will address the electronic structure of the 
FM states. Figure S] shows the total and atom-projected 
DOSs, which can be figured out well using the concept 



of RuOg cluster orbitals.pzj The Ru three t2g states, hy- 
bridized with the O six states, produce a bonding 
manifold at -2 to 1 eV, an antibonding manifold at -6 to 
-4 eV, and a nonbonding manifold at -4 to -2 eV, as can 
be visible in the lower inset of Fig. |4l The antibonding 
and bonding states of Ru two eg and O three states lie 
in the regimes of -8 to -6 eV and to 6 eV, respectively. 
The exchange splitting of the t2g manifold is about 0.5 
eV. In the FM state, the total moment of 1.01 /is is de- 
composed into 0.6 fiB for Ru and 0.28 for 30. (The 
remnant is in the interstitial regime.) 

The enlarged FM band structure of LSDA is given in 
the top panel of Fig. [Sj In this regime, the spin up and 
down structures are nearly identical, except for a differ- 
ence in on-site energy of roughly 0.5 eV due to exchange 
splitting. In the spin up, in addition to a flat band at 
-0.2 eV along the T-X line, which commonly appears in 
conventional perovskites, two saddle points at the X and 
M points occur, resulting in vHSs just below Ep. These 
features may cause the unusual electrical properties ob- 
served experimentally. [13] 

The band structure of the t2g and p^^ cluster in the 
range of -8 eV to 1 eV was reproduced well using the 
tight-binding Wannier function technique. We obtained 
two important parameters: pdir hopping t7r=l-23 eV and 
direct oxygen-oxygen hopping t'^—0.l7 eV. Compared 
with the tight-binding parameters of the nonmagnetic 
CaRuOs, obtained by Mazin and Singh, Q is about 
10% smaller, but t'^ is about half as large. This is consis- 
tent with the 20% smaller bandwidth of the antibonding 
t2g manifold in BaRuOa, implying stronger magnetic in- 
stability. 

Effects of spin-orbit coupling. SOC affects the Fermi 
surface topology in an isovalent Sr^Ru9O7.[30t though its 
effects are often negligible in Ad systems. To consider the 
effects of SOC, we conducted LSDA+SOC calculations. 
As shown in the upper inset of Fig. 21 the changes in 
the DOS are relatively small, leading to a tiny orbital 
moment of -0.013 fiB on the Ru ion. However, the mixing 
of states with other spin channel reduces the net moment 
by 10 %, bring it closer to the experimentally observed 
value of 0.8 fin- The net moment is 0.91 ^b, with 
local moments of 0.55 hb for Ru and 0.26 fiB for 30. 

An additional remarkable feature is that vHSs move 
toward Ep in both spin channels. In the spin up chan- 
nel, vHS at -0.1 eV splits into two peaks, whereas two 
vHSs merge and move to 0.1 eV in the spin down chan- 
nel. These features reflect changes in that band struc- 
ture at the X and M points, as demonstrated in the 
enlarged band structure in the bottom panel of Fig. [5] 
The strength of SOC is about 0.1 eV at these points near 
Ep, but negligible elsewhere. SOC alters the Fermi sur- 
face (FS) topology at the AT-points, in particular, in the 
intersecting pipe-like and cube- like FSs (not shown here) , 
which are often observed in conventional perovskites. [3l| 
The open pipe-like surface becomes closed, whereas the 
cube-like closed one becomes open. 
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IV. DISCUSSION AND SUMMARY 

The Fermi liquid behavior experimentally observed in 

the low T regime may suggest the presence of correlation 
effects. Considering that this is a metallic id system, 
however, the correlation effects should be weak. In our 
preliminary calculations using the LDA+U method, ap- 
plying U to the Ru ion shifts the Cg manifold up and the 
t2g manifold down, so only the intersecting pipe-like FS 
disappears. To evaluate the importance of the correlation 
effects, detailed experimental observations are necessary. 

In summary, we investigated the magnetic proper- 
ties of the ferromagnetic cubic perovskite BaRuOa, us- 
ing various first principles calculations: FSM, LSDA, 
LSDA+SOC, and the Wannier function technique. In 
the nonmagnetic DOS, the cubic phase has a sharp peak 
pinning at Ep, yielding IN{0)Kil.2, beyond the Stoner 
instability. In the LSDA calculations, a few vHSs appear 
near Ep iti both spin channels, resulting in a strong mag- 



netovolume effect of AM/Aaw 4.3 fJ^s/^, significantly 
higher than that of the Invar system ^Fe2. The effects 
of SOC are substantial near Ep around the X and M 
points and affects the topology of some of FSs, though 
the strength of SOC is negligible in most regimes. At 
the optimized volume, the effects of SOC reduce the net 
moment by 10% from the LSDA value, and it becomes 
very close to the experimentally observed moment. Our 
results show that the Stoner instability produces ferro- 
magnetism, and the combined effects of the magnetovol- 
ume and SOC determine the magnetic moment. Further- 
more, the unusually strong magnetovolume effects may 
imply atypical transport properties in this FM BaRuOs, 
suggesting further questions for experimental research. 
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